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Abstract. In the MSSM with complex parameters, loop corrections to the decay of a stop into a bottom 
quark and a chargino can lead to a CP violating decay rate asymmetry. We calculate this asymmetry at 
full one-loop level and perform a detailed numerical study, analyzing the dependence on the parameters 
and complex phases involved. If the stop can decay into a gluino, the self-energy and the vertex correction 
dominate due to the strong coupling. It is shown that the vertex contribution is always suppressed. We 
therefore give a simple approximate formula for the asymmetry. We account for the constraints on the 
parameters coming from several experimental limits. Asymmetries up to 25 percent are obtained. We also 
comment on the feasibility of measuring this asymmetry at the LHC. 



1 Introduction 

If the Minimal Supersymmetric Standard Model (MSSM) 
is realized in nature, LHC will produce squarks and gluinos 
copiously. However, even if supersymmetry is discovered, 
it will be still a long way to determine the parameters 
of the underlying model. In the general MSSM, the U(l), 
SU(2) and SU(3) gaugino mass parameters Mi, M 2 , and 
M3, the higgsino mass parameter fj,, and the trilinear cou- 
plings A f (corresponding to a fermion f ) may be complex, 



i.e. A f = \A f \ 



M 3 = \M 3 \e lv s. As usual, we take M 2 



positive and real by field redefinition pQ. The experimen- 
tal upper bounds on the electric dipole moment (EDM) 
of the electron, muon, neutron and several atoms severely 
constrain the phase of \i. In general, the phases of Mi, 
M3 and At b arc much weaker constrained due to possible 
cancelations [2l|51|il[51[rJl[71[51|gi[TfJl[TT|[T^ . Therefore, we use 
real values only for fj, and do not restrict the remaining 
phases. (For a recent discussion on EDMs see [13].) 

Complex MSSM parameters can lead to direct CP vi- 
olation (CPV), see the summary in [14]. One example is 
the CP violating rate asymmetry, which is a loop induced 
effect. CP violating asymmetries for the production and 
decays of the charged Higgs ff ± [IMIMnillHllISlISM^ 
|2"51T2"4] and for the decays ->■ W ± \° [25] were already 
studied in detail. Studies of measuring direct CP viola- 
tion in stop cascade decays at the LHC based on T-odd 
asymmetries built from triple products were done in [261 

mm- 

In the following, we study the CP violating decay rate 
asymmetry of the decays U — > b xl an d i* — > b x^ at mn 
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one-loop level [2"M50] . If the channel U — > g t is kinemati- 
cally open, the ti— t 2 self-energy and the vertex graph with 
g exchange are expected to dominate because of the strong 
coupling. But we show explicitly that the vertex contribu- 
tion is suppressed. All other contributions are numerically 
always smaller than ~0.5% which also means that the de- 
pendence on the phase of Mi is negligible. We thus give a 
short analytic formula for the decay rate asymmetry 5 CP 
which approximates the total one-loop result within 5% 
in the range above the threshold of the ti — > gt decay. 

In order to get a large decay rate asymmetry, not only 
the channel into g must be open, but also large phases 
or phase combinations of A t and M3 are necessary. In 
addition, the stops must be rather degenerate but with a 
strong mixing. The dependence on ipA b is weak because it 
only enters the vertex corrections. 



CP 



2 Decay Rate Asymmetry 6 



We define the CP violating decay rate asymmetry of the 
decays ti — > b xl an( i t* — > b x^ as 



5 CP = 



r+(u^bxt)-r-(t 



r+(u^bxi) + r-(t*^bx^)' 

The one-loop decay widths can be written as 

r± cx l-<cJ 2 + 2Re(^(X± co )tx i ± oop ) 

S S 

with the matrix elements given by 



(i) 



(2) 



Ml 



iu(k 1 )(B*P R + B^P L )v{~k 2 ) , 
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X treo = iu(k 2 )(B R P R + B^P^vi-h) , 
^loop- iu(h)(6B*P R + 5B L + P L )v(-k 2 ) , 
^loop = iu(k 2 )(6B*P R + 6BtP L )v(-h) 



(3) 



with the definitions 



^CP 



oc Re(6 g g 1 g 2 )Re(P&Ve) , 
oc — 2Im(6 go.9i.92)Im(PaVe) . 



(11) 



with P R x = (1 ± 7 5 )/2. The tree-level couplings B R ' L = 



B 



R,L 



B 



R,L _ r>R,L* 



B RL * are defined in Appendix iBl The form 
factors SB RL are calculated in Section [3J The form fac- 
tors 5B R ' L can be easily obtained by conjugating all the 
couplings involved. 

Since there is no CP violation at tree level, S CP is a UV 
convergent quantity which means no renormalization is 
necessary. Furthermore, we can write |A^^ CC | 2 as |.M trcc | 2 . 
Assuming that the one-loop contribution is small com- 
pared to the tree level, we use the approximation 



In order to obtain a non-zero 6 , not only the couplings 
but also the PaVe's must be complex. For that at least a 
second decay channel must be kinematically open, i.e. a 
particular one-loop diagram only contributes to the asym- 
metry if the corresponding two-body decay is kinemati- 
cally open. The asymmetry S CP = A^ p — A_ then be- 
comes 



S CP = [a{C rr + Ch L P ) - 2m b m^(C RL P + C£*)) 



/ElMrc 



(12) 



with 



S CP = 



r + -r- 

2-M;ree 



= A CP -A 



CP 



a cp _ M^s(Mt OD VMt oop ) 



(4) 



A(\B H \' 2 + \B L \ 2 ) - Am b m- + Re(B R *B L ^ 



and 



(5) 



Neglecting the bottom mass in Eqs. (|5I8I) the general for- 
mula of the decay rate asymmetry for a specific one-loop 
contribution simplifies to 



s cp^ 



1 



\B R f 



\B L \ 2 



FLR 
CP 



a 



LL\ 
CP) 



(13) 



Furthermore, we point out that in the asymmetry S CP 
possible rescattering effects (which are CP conserving) 
cancel each other and therefore drop out. 



Re(E(X± cc )tx i ± oop ) = A Re(B R 8B R + B%6BQ 



-2m b m i+ Re(B R SB^ + B^5B%) (6) 



using A = (m 2 — m 2 — m? + ). By defining combined cou- 

1 



pling matrices 



C l l = B^SBi 



(7) 



with i, j = R, L we obtain 
Re(£(M± ree )tjW± op ) = A (Re(Cr ) + Re(C^)) 

S 

-2m b m~+ (Re(C£ x ) + Re(C^)) .(8) 

These coupling matrices can be generally expressed by 
C± oc b± x (<?o<?i<?2)± x PaVe where b± = B^. is the cou- 
pling at tree level, (<7offi52)+ are the couplings of the three 
vertices and (<?o<7i<?2)- are the conjugated couplings. PaVe 
stands for the Passarino-Veltman-Integrals. Omitting the 
indices we can write 

Re(Ci') cx Re((6 ffo5lff2 ) ± x PaVe) = 
Re(6 50 3i32)Re(PaVe) t Im(fo. 90 5ig 2 )Ini(PaVe) (9) 

which leads us to the decomposition into CP invariant and 
CP violating parts 



(10) 



3 CP Violating Contributions 

In general, 47 one-loop diagrams can contribute. If the 
channel t j —¥ t g is kinematically open, the self-energy and 
the vertex graph (see Figure[T|) with gluino exchange dom- 
inate due to the strong coupling. The form factors 5B± 
are defined in Eq. ([3]). In the following, we only give the re 
suits for the form factor SB R since 5B+ 
L. The form factor for the self-energy process is 



SB R with R <-> 



SB' 



2 Cf m~ q m t 



o 

ru- 
t-i 



r>R 



(47r) 2 (m? 
■(G R *Gf + Gf*G R ^jB {ml 



2 2\ 

,m- g ,m t ) 



(14) 



with j 7^ i and Cf = 4/3. The form factor for the ver- 
tex correction reads (defining g^^ 1 = Gf*Gj*Al* with 
a,p,i = R,L) 



5B R 



C f 



2 



^ H L 

/ LRR i LLL 

\9ijk rn b m^++g ijk m g m t 
■9% RL m b m t + g$ R m-+m- g )C 



;i i (B (m?,m|,m?) 



RLL 



+ 771? Cq 



+(g!jk L m b + g RRL m g + 9^k R m it + g^ L m t )ra b C x 
+(gfjk R mb + 9it R m- g + g R i L m^+ + gf^ R m t )m^C 2 



(15) 
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Fig. 1. Feynman graphs with g ex- 
change contributing to S CP of U — > 

bjct — 1) 2; j 7^ i for the loop 
graph and j — 1,2 for the vertex 
correction graph) 



R,L 



B^ L , and G^' L are given 



The coupling matrices A k - 

in Appendix [Bj We use the one-loop integrals B , C , 
C\ , and C2 according to the definition of [3T] in the con- 
vention of [35]. The argument set for the C-functions is 
(to 2 ,to? ,to 2 + ,to? ,to?,to 2 ). 

In order to obtain 5 CP we insert 5B^ L (Eq. ([T4"f or 
Eq. (H5])) into Cf (Eq. ©). Then we calculate C^ p from 
Eq. (UH) and finally we get S CP from Eq. or Eq. (JT3J. 

We found that numerically only the self-energy graph 
is important. To understand the strong suppression of the 
gluino vertex graph in comparison to the gluino self-energy 
loop, one has to consider the possible combinations of the 
couplings in these graphs. 5 CP always contains a product 
of four squark rotation matrix elements. For the vertex 
graph they take the form R t R t *R b R b * , for the self-energy 
loop R t R t *R t R t *. Setting the indices of the external par- 
ticles to be t\ and xt i an d using the relations found in 
Eq. ([2"7} . we can rewrite these terms in terms of MSSM 
input parameters. 

Taking the input parameters given in Section [¥] we have 
mjj ~ m l 2 i m b 1 ~ m b 2 anc ^ gaugino like. Since we 
have mblm(Ab) -C m t lm(A t ) and assuming ip^ ~ 0, the 
remaining relevant term in S CP of the vertex graph is 
(A CP = g 2 C F m t lm(A t )/(ml - m? )) 



oc A CP g 2 



m xt 



2 

to- 

b 2 



2 

to- 

01 



(K 



m? )Im(C 2 (l)) 



v bl 



m? )Im(C 2 (2))) , 



(16) 



where A* 



is a trilinear breaking parameter, to? is the 



to-+ <C TOg and nearly degenerate sbottom masses. 

In the case of a xt which is higgsino like, the relevant term 
in SvJrtcx (which is now proportional to the Yukawa cou- 
pling \h t \ 2 instead of g 2 ) has a numerator which is again 
very small due to ~ and nib <C m t . Comparing 
this to the relevant term in 6^ (g 2 is simply replaced 
by — I /it 1 2 for hb <C h t ) one can see the same suppression 
mechanism at work. 

The suppression of the vertex correction is thus a general 
feature, even if xt becomes a mixed state. 

As a result for our scenario we can give an approxima- 
tive formula for S CP for the remaining leading self-energy 
contribution valid for r~ otal /2 <C \m*. — m?.l. Inserting 

Eq. (jnjl into Eq. (0 in order to calculate Eq. (JTTJ and 
finally Eq. UTS) we get 



5 CP 

approx 



1 



c f 



m~ g m t 



4n 2 (\B P \ 2 



\Bg\ 2 )(m 2 
■lm(b g g 1 g 2 )lm(B ) 

Im(6 5 o5iff 2 ) = Tm[(BgB% + B^B^) 



to? ) 



(18) 



Im(B ) 



using X(x,y,z) 



ttA 1 /2( 



2 9 Ox 



= a; 2 + y 2 + z 2 

jtCP 



0(™ t -. 



2xy — 2x2 



TO t ) 



(19) 



step function 0. 5H,p rox is a good approximation of S^f (all 



(1,1) element of the sbottom mass matrix, and C 2 (j) is contributions) above the threshold of the t 
the Passarino-Veltman integral with m^ . In the limit 



2yz and the 
Lof^f (al 
g t decay. 



mg — > to^ 2 (and thus 62(1) — >• C*2(2)) <S^,f tcx vanishes. 
But even if the sbottom masses would not be degenerated 
(and thus yielding a higher numerator) the denominator 
always compensates this effect. The relevant term of the 
self-energy loop is 



4 Numerical Results 



xCP 

-'self 



A CP g 2 



2 

mj 

12 



2 

mj 



■Im(B 



(17) 



Comparing S^' p tcx with 5^ we can see that the suppres- 
sion of the gluino vertex correction is due to mb <C m t , 



We present numerical results for the decay rate asymme- 
try 8 CP as well as the tree-level branching ratio (BR) of 
the process t\ — > bxt- The 47 one- loop contributions to 
S CP were calculated by using FeynArts [33] . Further- 
more, the gluino graphs and a few other ones were cal- 
culated independently and also cross checked numerically. 



For the Im(So) relation see Eq. (79) in [21] 
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The strong coupling a s is taken running in the DR scheme 
at the scale toj . In the calculation of 5 CP we also take 
the Yukawa couplings (h t , hb) running as given in the Ap- 
pendix A of [15] , 

We calculated the EDMs up to leading two-loop order 
with CPsuperH [33] to check that our parameter points 
are consistent with the constraints coming from the EDM 
of the electron, muon, neutron (all [35] V and mercury [36] . 
We can fulfill these constraints since we take n real, choose 
only the third generation breaking parameters A t .b, T to be 
complex, and because we can always choose the squark 
SUSY breaking parameters of the first and second gener- 
ation appropriately. We get the right amount of the cold 
dark matter relic density [37J with the Xi LSP annihilating 
mainly into r+ t~ (micrOMEGAs [3EUM]) by varying 

£ so that m^o ~ m,f x . Furthermore, the constraints 
coming from B — > X s 7, B s — > fi + fj," and Bd — > t + t~ 
(all [40]) as well as the Higgs mass limit [35] are fulfilled. 

For the numerical analysis, we fix for the third gener- 
ation Mq = M = Mjj, M L = M& \At\ = \A b \ = \A T \, 
\Mi\ = M2/2, and for the complex phases fA t = fA b = 
fA T - We start from the following MSSM reference sce- 
nario: Mq = 650 GcV, M L = 120 GeV, \A t \ = 190 GeV, 

<f At = tt/4, = <PMi =^g = 0, M 2 = 150 GeV, 
\fi\ = 830 GeV, tan/3 = 5, and M A o = 1000 GeV. These 

412 GeV, m~ ti = 
0+ 



parameters give m~+ 



146 GeV, rrig 



653 GeV, m t ~ 2 = 688 GeV, r|° tal = 18 GeV, xi is gaug- 

ino like, and t\ and ti have a low mass splitting but large 
mixing. 

In Fig. [3 and Fig. [3] we show 8 CP and BR as a function 
of mj for f A t = f i T"- Higher values of fA t result in a 
less degenerate stop mass splitting which reduces the en- 
hancement of S CP coming from the stop propagator (see 
Eq. (|18pl. For the parameter Mq is varied from 500 

to 1000 GeV. One can see the threshold of the t% gt 
decay at ~583 GeV, after which the gluino contributions 
account for up to ~ 98% of 6 CP . However, if this decay 
channel opens, the BR drops quickly. This feature can 



ti -> bxt 




500 600 700 800 900 1000 

m~ h [GeV] 

Fig. 2. 8 CP as a function of (Mq varied) for various values 
of f At 




500 600 700 800 900 1000 

m h [GeV] 

Fig. 3. BR as a function of (Mq varied) for various values 

Of f At 



is large but BR is small. The other decay channels are 
h -> txtbxt- 





80 : 




60 - 




40 - 




ft? 




cq 


20 - 







ti -^bxi, ti 
i i i i i i i i j 



gt 



600 



1000 



be seen in detail in Fig. [U If t% —> g t is possible, 5 



CP 



700 800 900 

m k [GeV] 

Fig. 4. Comparison of BR(ii — > &xf ) an d BR(ii gt) as a 
function of {Mq varied) 



We also studied the dependence on the gluino phase <p§ 
as a second source of CP violation, which is, however, in 
conflict with the current EDM limit of mercury by a factor 
of two. 

Figure [5] shows the contributions from the self-energy and 
vertex graphs with gluino exchange (see Fig. [1]) as a func- 
tion of mi . As already anticipated in Section[3J the gluino 
self-energy loop dominates. 

In Fig. [6] we show the ratios of 5 between the ap- 
proximated formula for the gluino self-energy loop <5^pp rox 
(Eq. (fT5)l) and all one-loop contributions S^iC > as wen as 
between both gluino contributions 8^ p and all contribu- 
tions. Above the threshold of t\ — > gt both gluino pro- 
cesses account for ~ 97% of all processes. One can see 
that <5fpp rox is indeed a good approximation of S^f . 
We also studied the dependence of S CP on tan/3 and \A t \. 
The main effect comes from the off-diagonal element in the 
stop mass matrix a t = A* t — /x/tan/3. For \A t \ ~ 190 GeV 
and tan ft ~ 6 the asymmetry has its maximum up to 
25% because a t becomes minimal. In this case one has 
rather degenerate stop masses which enhance the gluino 
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600 



700 



800 



900 



1000 



m ix [GeV] 



Fig. 5. Contribution of the gluino self-energy loop (1) and the 



vertex correction (2) to 8 CP 

ii - 



as a function of (Mq varied) 




600 700 800 900 

m h [GeV] 



1000 



Fig. 

A- a -. 



6. Ratios of A 
sCP i x CP 
3 ~s 



8 CP 



approx 

8^ P /8^T as a function of mj 1 (Mq varied) 



/8%F, see Eq. (TH}, and 



self-energy contribution due to the propagator oc 1/ (m~ — 
m^). For larger values of \A t \ and tan/3 the asymmetry 
decreases and begins to be in conflict with the EDM limit 
of mercury 

The effect on the mass splitting of t\ and t% is shown in 
more detail in Fig. [7] and Fig. |5J We fix mj = 650 GeV 
and vary by changing the parameters Mq, Mfy. There 

are two possibilities, toll §s mail- In our scenario xt * s 
gaugino-like which couples dominantly to t^- For toll < 
ttirr we have t\ ~ ti, and t% ~ Ir. Hence BR(i\ — > bxt) 
is large, but S CP is small because the coupling of the in- 
ternal ?2 to bxt i s suppressed (see left graph of Figured]). 
For toll > mRR one has the opposite behavior. When Xi 
is higgsino-like the whole situation is reversed. We there- 
fore see that for the cases where S is large the BR is 
always small and vice versa, unless the masses of t\ and 
t-i are rather degenerate. Note that the masses cannot be 
arbitrarily degenerate, since otherwise the off-diagonal el- 
ements in the stop mass matrix (and thus the complex A t 
parameter as the main source of CP violation) would need 
to vanish. 

Furthermore, it is crucial that the enhancement of S CP 
due to a small stop mass difference is not a resonance en- 
hancement, i.e. the mass difference mj 2 — mj must be 
sufficiently large compared to the FWHM r~ otal /2 of the 



resonance. In our scenario this is always the case, since 
rtotal i 2 = 9 GcV and A m = 35 GeV a t the reference 

point and r^ otal /2 always remains much smaller for vari- 
ous mass differences. 




700 750 800 850 900 950 1000 

m h [GeV] 

Fig. 7. Effect on mass splitting of ti and I2 for S CP 
function of rrij 2 (Mq , Mq varied) for 



650 GeV 



'XT 



30 
25 
_20 
"15 

5 




m LL < m RR 



m LL > m RR 



Fig. 8. 

tion of 



700 750 800 850 900 950 1000 

m h [GeV] 

Effect on mass splitting of t\ and ti for BR as a func- 
m t - 2 {Mq,Mjj varied) for m il = 650 GeV 



For completeness, we also examined the process ti — > b xt ■ 
Because of the similar masses and large mixing of the 
stops, the resulting plots are alike. We obtain S CP = 19% 
and BR = 17% at the reference point. 
We also investigated the influence on the Yukawa cou- 
plings h t and /if, taken to be running. In our scenario, the 
difference of the asymmetry S CP taken with running and 
non-running Yukawa couplings is negligible. 
The theoretical uncertainty of the asymmetry S CP is esti- 
mated ~20%, due to higher order corrections [HH22]. 
For a measurement of the asymmetry 5 (Eq. (fT])) at 
LHC, one has to consider the process 



pp -> hh + X -> (bxt){bxl) + X 
-> (bW+xl)(bW-x°i) + X 
-> 2 b-jets + W+ + W-+ $ T + X 



(20) 
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(where X only contains the beam jets) with one W + de- 
caying hadronically and the other one leptonically to get 
information on the charge of the W, and in turn of the 
chargino and the stop. We assume that the masses of t\, 
?2 1 xf and X? are known so that one has enough kinemat- 
ical constraints to single out the respective decays [33] . In 
such a way it should also be possible to separate off the 
production and decay of t\ from those of t^. For yfs — 

14 TeV, mjj = 610 GeV the cross section of pp — > t\t\ 
is (7 = 200 fb at NLO according to Prospino [33]. As- 
suming C = 300 ftT 1 and BR^ -S- b X t) = 0.2 we esti- 
mate a purely statistical relative error of the asymmetry 
AS CP /6 CP = ±0.077. The largest background is of course 
pair production of top quarks [35] . For a realistic estimate 
of the measurability a Monte Carlo study would be neces- 
sary, which is, however, beyond the scope of this article. 
For the same decay chain of t\ such a study was done for 



an e T e 



collider 



For the stops (/ = t) we simply insert the corresponding 
values 7 t 3L = 1/2, e t = 2/3, Mq, Mq, m t and A t . Anal- 
ogously, for the sbottoms (/ = b) we have I^ L = —1/2, 
ei, = —1/3, Mq, Mjj, mi, and At,. 

M. 2 j is diagonalized by the rotation matrix BJ such that 



RfM 2 j{Rfy = diag(mj-,mj 3 ) and (J 
have 



Rf{ 



In' 



We 



uf - ( r il r ir 1 _ 
\R f R f 



e IV /sinl 



-e r f smdj cos 6^ 



(25) 

Using the unitarity property of the sfermion rotation ma- 
trices and the diagonalization equation for the sfermion 
mass matrix, 



M ) = m K mf : f 



(Rfy [o 1 ™}J Rf > (26) 



5 Conclusions 

In the MSSM with complex parameters, loop corrections 
to the ii — !• fex/t decay can lead to a CP violating decay 
rate asymmetry S CP . We studied this asymmetry at full 
one-loop level, analyzing the dependence on the parame- 
ters and phases. Below the threshold of the ti — > gt decay, 
S CP is < 1%. If this channel is open, a 5 CP up to 25% 
is possible (mainly due to the gluino contribution in the 
self-energy loop), if the stop particles have a small mass 
splitting together with large mixing and the chargino is 
wino like. 
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A Masses and Mixing Matrices 

The sfermion mass matrix in the basis (/l, fn) with / 
t,b,f,... is 



Mi =i 2 
/ y mf a f mj R 

with the following entries 
.- _ u- 

,2 „„„oo/t3L „ „;„2 



21) 



m h = M {Q;L}+ m ? 



+mz cos 2(3(1 j - e f sin 2 W ) 

2 ti /r2 | /m 2 I „_ 2 n n • 2 . 



(22) 



mj R = M0.f).^y +mj + m z cos 2/3 e/ skr 9 W , (23) 



J m u (>l* — n cot /?)... up- type sfermion , 
' 1 m <i{A* d — //tan/3) ... down-type sfermion . 



one can derive the following relations (we define A$ = 
mml-ml)): 

R f 22 R f 2 l = -R{ 2 R{1 = m f a f Af, 

R f 21 R f 2 *2 = -RiMl = m f a}Af, 

RIxRII = R{2 R il = (m% - mjJAf , 

R( 2 R{* = R^Ril = -{m\ ~ m 2 h )A* . (27) 

The chargino mass matrix in the basis (— iX + ,ipj l2 ) is 



M c = 



Mi \/2 mw sin /3 



(28) 



(24) 



mw cos /3 fi 
It is diagonalized by the two unitary matrices U and V 

U*M C V^ = diag(m-± ,m~±), (29) 
where m -± arc the masses of the physical chargino states. 

Xl,2 ^ J b 



B Interaction Lagrangian 

In this section we give the parts of the interaction La- 
grangian that we need for the calculation of the lead- 
ing contributions. The chargino-squark-quark interaction 
is described by 

*W = t( A kiPR + A L kl P L )xt k 

+b[B p P R + B^P L )xt c U 

+^A L k *P R + A p *P L )tb* 

+W(B%P R + B?*P L )b t* . (30) 

The couplings are 
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V2 \mw cos p 



gmt 



* rib* 



-VCoR 



k2 — fK ■ n k2 1 Hl > 

V 2n 



2mw sin j3 

B& = KR%v k 2 - gRtv kl 

= 4= ( m \ V k2 R% - V2V kl Rll) , 
y/2\m W sm0 J 

B L = h bRTi U k2 = 7k 9 b o U k2 R n > ( 31 ) 

V Zmw cos p 

where we used the relations m t — (h t ^sin/3)/v / 2, mb — 
(hb ucos0)/\2 and raw = to convert the coeffi- 

cients. 

The gluino-squark-quark interaction is 

+F(G%£ P R + G^: P L )g a q? , (32) 

where we used u, v as the color index (u, v = 1, 2, 3), i as 
the mass index (i = 1, 2) and a as the gluon/gluino index 
(a = 1, . . . , 8). The couplings are 

= T« Gf with Gf = V2.g s i?f 2 , 
G^ v =T^Gf withGf = -V2g s Rl 1 e- i ^ (33) 

with T" v as the generator of the SU(3)c group, g s as the 
coupling constant of the strong interaction and ipg as the 
gluino mass phas^E 
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